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Q ' Abstract 
(D ! 

Wc consider nonlinear biasing models of galaxies with particular attention to a correlation between 
linear and quadratic biasing coefficients, hi and 62- We first derive perturbative expressions for hi and &2 
CO ' in halo and peak biasing models. Then we compute 

power spectra and bispectra of dark matter particles and halos using N-body simulation data and of 
volume-limited subsamples of Sloan Digital Sky Survey (SDSS) galaxies, and determine their hi and 62. 
We find that the values of those coefficients at linear regimes (fc < 0.2ft,Mpc~^) are fairly insensitive to the 
' redshift-space distortion and the survey volume shape. The resulting normalized amplitudes of bispectra, 

\ Q, for equilateral triangles, are insensitive to the values of hi implying that ^2 indeed correlates with hi. The 

. present results explain the previous finding of Kayo et al. (2004) for the hierarchical relation of three-point 

0^ ' correlation functions of SDSS galaxies. While the relations between hi and &2 are quantitatively different 

, for specific biasing models, their approximately similar correlations indicate a fairly generic outcome of the 

biasing due to the gravity in primordial Gaussian density fields. 
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?H ] and understanding is crucial in advancing our knowledge of the evolution of the universe. 

A reasonable approximation often adopted is a local linear biasing model, which assumes a simple scale-independent 
relation between the density contrast fields of galaxy and mass: 

5g{^,z)^h{z)6{^,z). (1) 

In the above, the density contrast of the mass is defined as 

5(x,.).^i^^li|-^, (2) 

where the over-bar indicates the mean over the entire universe. The number density contrast of galaxy density field, 
5g, is defined similarly. Because this model contains only a single time-dependent parameter h{z), it is not surprising 
that the model cannot describe the observed features of galaxy clustering accurately. Furthermore, the three-point 
correlation function of SDSS galaxies (Kayo et al. 2004) indicates the importance of the higher-order correction to 
equation (1); while the galaxy two-point correlation functions are well represented by the linear biasing model (e.g., 
Zehavi et al. 2005), the normalized amplitudes Q of the corresponding three-point functions for equilateral triangles 
do not show the expected scaling with respect to b(z): 

Q«i. (3) 
In reality, however, Q for equilateral triangles calculated from SDSS galaxies in redshift space proves to be almost 



1. Introduction 

One of the major uncertainties in precise cosmology is the galaxy biasing relative to the underlying mass distribution. 
It hampers extracting the cosmological information from the large scale structure of the universe. In particular, the 
biasing is sensitive to the unknown physical conditions of galaxy formation. Thus its phenomenological parametrization 
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scale- independent, and follows the hierarchical relation approximately, Q = 0.5 ~ 1.0. Moreover its dependence on 
the morphology, color, and luminosity is not statistically significant. Given the robust morphological, color and 
luminosity dependences of the two-point correlation function, Kayo et al. (2004) argued that galaxy biasing is complex 
and requires a contrived relation between the linear biasing and its higher order correction terms. Another possibility 
is that the observed value is significantly contaminated by the redshift-space distortion and does not properly reflect 
the actual clustering information in real space. Indeed the previous N-body studies of the redshift-space distortion on 
three- and four-point correlation functions suggest that this is the case in nonlinear regimes (Suto 1993; Matsubara 
& Suto 1994; Suto & Matsubara 1994). The interpretation is further complicated by the fact that the higher-order 
statistics is sensitive to relatively rare large-scale structures in particular samples (Nichol et al. 2006). 

The main purpose of the present paper is to provide a physical explanation for the approximate hierarchical relation 
for Q of SDSS galaxy clustering on the basis of perturbation analysis and numerical simulations. This is a step toward 
the understanding of the nonlinear nature of biasing. Specifically we employ a local nonlinear biasing model (Fry & 
Gaztanaga 1993): 

5,(x,z) = f;^[<5(x,z)r. (4) 

We consider density peaks, dark matter halos, simulated halos, and the SDSS galaxies as specific examples for the 
density fields 5b in the left hand side. 

The outline of the paper is as follows; section 2 briefly summarizes the halo and peak biasing models as analytically 
tractable examples. We compute the linear and quadratic biasing coefficients perturbatively and find that these 
models roughly agree with the observed hierarchical relation. To be more realistic, we analyze simulated halo catalogs 
and SDSS galaxies in section 3, and find that the pcrturbative result is valid even if we take account of a variety of 
selection effects and particularly redshift-space distortion. The quantitative agreement between the two indicates a 
generic correlation among the biasing coefficients in the SDSS galaxies. Finally, section 4 is devoted to the summary 
and conclusion of the paper. 



2. Perturbative predictions in halo and peak biasing models 



2.1. Basic statistics and perturbation biasing models 

Throughout the present analysis, we work in the Fourier space. For definiteness, we adopt the definition of the 
Fourier transform of an arbitrary function ^(x) as 

'^^^ " / (^^We-"^-. (5) 
Then the power spectrum and the bispcctrum of the mass density field are defined as 

(5(ki )5(k2 )) = P(fci )5d (ki + k2) , (6) 

(<5(ki)<5(k2)<5(k3))=S(fci,fc2,fc3)<5l3(ki+k2+k3), (7) 

where 5d denotes the Dirac delta and we assume the universe is isotropic and homogeneous. We introduce the 
normalized amplitude of the mass bispectrum as 

n (1, h lc\^ B{ki,k2,k3) 

"^"'^ P{k,)P{k2) + P{k2)P{ks) + P{ks)P{k,)- 

Note that strictly speaking, the above statistic is different from that defined in configuration space^. Since both are 
expected to be nearly identical in a weakly nonlinear regime, we focus on the Fourier space analysis throughout this 
paper (see also Verde et al. 2002 and Hikage et al. 2005). The analysis in configuration space is now in progress. 

The above statistics can be generalized in a straightforward manner to any biased field. If one keeps the linear and 
quadratic terms in equation (4), one obtains 

Pb{k)^blP{k), (9) 



Qb{ki,k2,k3) ^ 

Ol 



Qm(fcl,fe,fc3) + 7^ 
Ol 



(10) 



Therefore, Qf, for the biased field is directly dependent on the quadratic biasing coefficient 62 unlike in the case of 
the power spectrum. As we show below, this term plays an important role to approximately satisfy the hierarchical 



In this paper, we use the term "real space" in order to imply the analysis without redshift-space distortion in k-space. Thus we reserve 
the term "configuration space" in order to distinguish from k-space. 
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relation. In order to understand the qualitative relations among the biasing coefficients, we first consider analytically 
workable biasing models, i.e., halo and peak biasing. 

2.2. Halo biasing 

We follow the formalism of Mo & White (1996) for halo biasing. We define (unconditional) mass function, nhaio('Ti,^) 
of halos with mass m at redshift z. We rewrite the unconditional mass function as 

^vf{v) — , v(m,z) = dsc[z)/a{m). (11) 

We denote by CT(m) the linearly extrapolated value of the rms of the initial density fluctuation field smoothed with a 
tophat filter of scale R = (3m/47rp)^/'^, where p is the mean comoving background density. The characteristic density 
contrast for spherical collapse is given as 
3('127r^2/3 

Ssciz)- ^ [l + 0.0123log,Mz)]/D{z), (12) 

where flmiz) is the matter density in units of the critical density at redshift z, and D{z) is the linear growth rate 
(Kitayama & Suto 1996). In the above equation and throughout the paper, we assume a spatially-flat model with non- 
vanishing cosmological constant. Specifically, we adopt the present value of mass density ilm = 0.3, the cosmological 
constant J^a = 0.7, the Hubble constant h = i7o/(100km s~^Mpc~^) = 0.7, and the amplitude of the density fluctuations 
smoothed over 8/i~^Mpc erg = 0.9. 

There exist two popular models for /{ly). One is based on a spherical collapse model (Press & Schcchtcr 1974): 

z./ps(z^) = 2Jj^exp(-;/V2). (13) 
The other is based on an ellipsoidal collapse model (Sheth & Tormen 1999; Sheth, Mo & Tormen 2001): 

z./st(^) = 2A{p) (^1 + j-^^ y^exp(-q^V2), (14) 

where 0.3, A{p) = [1 + 2-Pr(l/2 -p)/V^-i w 0.3222, and g« 0.75 (Cooray & Sheth 2002). 

Consider a comoving volume V and its total mass M at redshift zq. Let nhi^io(m, zi\AI,V, zo)dm denote the con- 
ditional number density of halos whose mass is between m and m -I- dm at zi . A reasonable approximation for 
n-haioiiTi, zi\M ,V, Zq) IS obtained by applying the extended Press-Schechter theory (e.g.. Bower 1991; Bond et al. 1991): 

m'^nh^ioim,zi\M,V,ZQ) dm dv^ 

= = , (15) 

P m 1^10 

Sscjzi) - Sojzo) 

= — , (16) 

^a^{m)-a^{M) 

where 5q{zq) is the linearly extrapolated mass density contrast at zq. If we define the actual density contrast 5 = 
M/{pV) — 1, 5q[zq) is written as 

^^g.^S (17) 

17 341 55805 

ao = 0; ai = 1; a, = --; a, = «4 = -y^; (18) 

where the above coefficients are obtained assuming the spherical collapse model (Bernardeau 1994). 
Now one can write down the corresponding conditional density contrast of halos: 

X ( ifl/fi/ \ ni,^io{m,zi\M,V,Zo){l + 5) 

Ks.io{m,zi\M,V,ZQ)^ 1. (19) 

"-halo(?Tl,Zl) 

Expanding equation (19) in terms of 5q in the limit of a{M) yields the biasing coefficients of halos (Mo, Jing & 
White 1997; Scoccimarro et al. 2001; Cooray & Sheth 2002): 

bi{m,z) = 1+61 + El, 

62(to,z) = 2(l + a2)(ei+£;i) + e2 + £^2, (20) 
where 
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Fig. 1. biasing coefficients in the three perturbative biasing models. Thick and thin curves show bi and 62/^1 at 
z = (soUd) and z = 2 (dashed). We use spherical halo model (left), ellipsoidal halo model {center), and peak model 
[right). The biasing coefficients are plotted against halo mass Mhaio (left and center), and against peak height 1/ (right). 



El 



2p 



'P 



2ei. 



(22) 



The above results for ellipsoidal collapse reduce to those for spherical collapse, if one sets p = and q=\. 
2.3. Peak biasing 

The biasing coefficients of the peak model can be obtained similarly as the halo model described in the above 
subsection. The conditional and unconditional number densities of peaks with peak height v = bja are derived in 
Bardeen et al. (1986). In this case, a denotes the rms value of density fluctuation smoothed over R = (iraj ^ixpYI"^ , 
and 8 is the density contrast of the peaks in the smoothed field. Substituting those formulae into the right-hand side 
of equation (19), one similarly obtains the biasing coefficients for the peak model (Mo, Jing & White 1997): 



6i(i/,z) = l + 



31 



<5sc(2:) 



62(1^,2) = 2(1 + 02)- 



.91 



2 1,0 I 
V - 1 + 2,gi H ^ 



(23) 
(24) 



(5sc(2) \8sc{z) ^ 

where the above functions, g\ and gi, are defined in equation (25) of Mo, Jing & White (1997). 

2. 4. Results of Analytic Models 

We consider the above three models, spherical halo, ellipsoidal halo, and peak, to explore the correlation between 
bi and &2 analytically. 

Figure 1 shows biasing coefficients, 61 (thick lines) and 62/^1 (thin lines) at z = (solid) and z = 2 (dashed). We adopt 
the cold dark matter transfer function of Bardeen et al. (1986) as the initial mass power spectrum (with ri,„ = 0.3, 
nA = 0.7, and the shape parameter T = 0,„/i = 0.21). In the three biasing models, 61 and 62/fei behave similarly as 
functions of Mhaio or i>. This implies a presence of a certain correlation between bi and &2/&1- To see this point more 
clearly, we plot 62/^1 in terms of 61 for each model (figure 2). All the three models exhibit very similar correlations 
between 61 and b2/bi. The differences among the three models become even smaller at higher redshifts; compare z ~0 
(thick lines) and z — 2 (thin lines). 

In analyzing simulation data below, the estimate of the biasing coefficients is made after averaging over a finite mass 
range. We calculate the mass-averaged values of bi and &2 for halo models at z = 0: 



Eji (^min ; ^max ) 



dmn]^a\a{m,z ~ 0)&„(m,z = 0) 



(n = l,2). 



(25) 



dmn\,s.\o{'m,z = 0) 

The difference between &„ and S„ is illustrated in figure 2, where spherical (open) and ellipsoidal (filled) halo models 
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Fig. 2. Correlation of bi and 62 /t'l from halo and peak biasing. Different line types are results of halo biasing of spher- 
ical collapse{solid), ellipsoidal collapse(dashed), and peak biasing (dotted), evaluated at redshifts 2 (thin) and (thick). 
Open (filled) symbols represent the mass-averaged values, Bi and B2/B1, defined in equation (25), assuming spheri- 
cal (ellipsoidal) halo model. The mass ranges correspond to L (square), and S (triangle) defined in subsection 3.1. 
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Fig. 3. Q for biased fields as a function of h\ in real space. We estimate Q using equation (10) on the basis of per- 
turbative predictions for peaks (dotted), ellipsoidal halos (dashed), and spherical halos (solid). We assume three different 
values for Qm'- Qm = 1.58, Qm = 0.83, and Qm = 4/7 from top to bottom. The first two values are computed from 
N-body results at A; = 0.4/iMpc~^, and k = 0.18/iMpc~^ (see section 3 below), and the last value corresponds to pertur- 
bation theory [eqs.(39) and (40)], where equilateral triangles (fci = k2 = ks = k) are assumed. left: z = 0, right: z = 2. 



[eqs.(13), (14)] are assumed. We plot the results for two mass ranges, corresponding to S (triangle) and L (square) 
halo subsamples (see subsection 3.1 for detail). 

Figure 3 plots Qb for biased fields as a function of bi in real space. The dotted, dashed and solid lines represent 
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Fig. 4. Same as figure 3, but for rcdshift space. We use equation (26) for equilateral triangles instead of equation (10). In 
this time, Q\, is calculated directly (not in terms of Qm), so we plot only one set corresponding to Qm = 4/7 in figure 3. 



the results of tree-level perturbation with respect to biasing, equation (10), for peak, ellipsoidal, and spherical halo 
models. For simplicity, we consider equilateral triangles and adopt three different values for Q^', Qm = 4/7 from the 
leading term of the gravitational nonlincarity [see eq.(39) below], and Qm = 0.83 and 1.58 from N-body results at 
k = 0.18/iMpc"^ and k = 0.40/iMpc"^ 

Figure 4 plots Qb for biased fields but in rcdshift space. The redshift distortion effect for Q is difficult to model 
accurately, and wc follow Scoccimarro, Couchman & Frieman (1999) and replace equation (10) by 

^ 5(2520 + 44IO7+ 1890/3 + 29407/3+378/32 + 4417/3^ + 9/3^ + 1470&i/3 + 88261/3^ - Uhf]"^) 

Qb,scq = noL /i c , i n ,3 , n a2\2 ' (^O) 



986i(15+10/3+3/?2)2 



"4 



4/7 SiA 

™ + 70 



1 



(27) 



for equilateral triangles, with 7 = 62/^1. 

These simple halo and peak biasing models imply that the nonlinearity of biasing weakens the dependence of 
Qhaio and Qpoak On 61, which qualitatively explains the results of Kayo et al. (2004) for equilateral triangles. To be 
more realistic, however, we consider in the next section nonlinear gravity, rcdshift-spacc distortion and finite survey 
boundary effect using N-body simulations and SDSS data. 



3. Comparison with N-body simulations and SDSS galaxies 



3.1. N-body simulations and SDSS galaxies 

N-body simulations that wc use employ 512^ dark matter particles in a cubic box of 300/i~^Mpc (comoving) on a 
side with the periodic boundary condition. The mass of each dark matter particle is 1.68 x lO^"ft.~^M0. As in section 
2, we use the cold dark matter transfer function of Bardeen et al. (1986) as the initial mass power spectrum (with 
r2,„ = 0.3, Q\ = 0.7, and F = ^mh = 0.21), and have generated three independent Gaussian realizations. The initial 
density field at z 36 is evolved up to z = using the P^M code of Jing & Suto (1998, 2002) with the gravitational 
softening length of e ~ 58/i~^kpc. We assume the amplitude of the density fluctuations smoothed over 8ft.~^Mpc 
as = 0.9. 

In order to test the dependence of Qb on actual galaxy properties, wc construct volume-limited subsamplcs of 
different colors from New York University Value-Added Galaxy Catalog (Blanton ct al. 2005) based on the SDSS 
Data Release 4 (Adelman- McCarthy et al. 2006). The angular selection function of the survey is written in terms 
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Table 1. Volume-limited subsamples of SDSS galaxies for 0.05 < z < 0.1. 



subsamplc 


luminosity 




color 






27rny^[/iMpc"^] 


bi 


blue 


-21.3 < M,. < - 


19.8 


g-r< 0.86 


33986 


0.00359 


0.962 


0.95 


red 


-21.3 < M, < - 


19.8 


g-r > 0.86 


34351 


0.00363 


0.966 


1.33 



Table 2. Subsamples of simulated halos. 



subsample 










27r(nh)i/3[/iMpc"'] 


ih) 


r-cubc S 


1.68 X 10^^ — 1.18 X 


10" 


58692 


0.00217 


0.814 


0.83±0.01 


r-cubc L 


> 1.18 X 10^3 




10881 


0.000403 


0.461 


1.26±0.02 


r-cubc LL 


> 6.72 X 10^3 




1484 


0.0000550 


0.239 


1.8±0.1 


s-cubc S 


1.68 X 10^2 — 1.18 X 


10" 


58692 


0.00217 


0.814 


0.88±0.01 


s-cubc L 


> 1.18 X 10^3 




10881 


0.000403 


0.461 


1.29±0.01 


s-cubc LL 


> 6.72 X 10^3 




1484 


0.0000550 


0.239 


1.8±0.1 


wedge S 


1.68 X 10^2 — 1.18 X 


10" 


21735 


0.00230 


0.829 


0.84±0.01 


wedge L 


> 1.18 X 10" 




3979 


0.000421 


0.471 


1.24±0.09 


wedge LL 


> 6.72 X 10" 




551 


0.0000582 


0.243 


1.8±0.2 



of spherical polygons (Hamilton & Tegmark 2004). Details of the SDSS can be found in the following papers: York 
et al. (2000) describe an overview, technical articles providing details include descriptions of the telescope design 
(Gunn et al. 2006), the photometric camera (Gunn et al. 1998), photometric analysis (Stoughton et al. 2002), the 
photometric system and calibration (Fukugita et al. 1996; Hogg et al. 2001; Ivezic et al. 2004; Smith et al. 2002; Tucker 
et al. 2006), the photometric pipeline (Lupton et al. 2001), astrometric calibration (Pier et al. 2003), selection of the 
galaxy spectroscopic samples (Eisenstein et al. 2001; Strauss et al. 2002), and spectroscopic tiling (Blanton et al. 2003). 

We only consider galaxies with rcdshifts 0.05 < z < 0.1 and magnitudes —21.3 < Mr < —19.8 (68337 galaxies in 
total, over ^ 4259dcg^). We divide those galaxies according to their colors so that each subsample contains roughly 
the same number: red subsample has g — r > 0.86 (34351 in total), and blue subsample has g — r < 0.86 (33986). In 
Table 1, Ng and Ug denote the total number of galaxies and the galaxy number density in each subsample, and the 

1/3 

wavenumber corresponding to the mean galaxy separation is 2TTng . The linear biasing coefficient, 6i, is evaluated 
at = 0.126/iMpc~^ assuming the same set of cosmological parameters used in N-body simulations. In particular the 
estimated value of bi is sensitive to as. The value of cts is still in controversy between different observations. For 
example WMAP3 gives as =0.742 ±0.051 (Spergel et al. 2006), while 2dFGRS gives 0.88+° (Gaztafiaga et al. 2005). 
Thus in principle our estimated values of bi would increase by about 20% if we had adopted WMAP3 result. 

In order to compare with analytical predictions for halo biasing and also with the SDSS subsamples, we identify 
dark matter halos using a fricnds-of-friends algorithm. Specifically we use the public code "FOF" ^, and adopt 
the linking length 0.164 in units of the mean particle separation. We construct nine halo subsamples from three 
realizations (Table 2). The indices, S, L, and LL, represent the three different halo mass ranges corresponding to 
1.68 X lO^^/i-iMg < M < 1.18 X lO"/i-iM0, 1.18 x IO^^H-^Mq < M, and 6.72 x lO^^h'^MQ < M. The mass ranges 
of the first two, S and L, are determined so that they have approximately the same values of bi for the SDSS blue 
and red galaxies, respectively. The last one, LL, is constructed so as to check the 62-61 correlation around 61 ~ 2, 
and discussed in figure 8 below. We further divide the three different mass-selected samples into three subsamples: 
r-cube (s-cube) measures positions of halos in real (redshift) space using the original simulation cube, (3006.~^Mpc)'^. 
In the latter, redshift-space distortion effect is taken into account by using z-component of the center-of-mass velocity 
of each halo. Wedge subsamples are constructed so as to have the same survey geometry as the SDSS sample. They 
measure the positions of halos in redshift space, using the line-of-sight component of the center-of-mass velocity unlike 
s-cube. In Table 2, (Nh) and (n^) denote the total number of halos and the halo number density in each subsample 
(averaged over three realizations). The wavenumber corresponding to the mean halo separation is 27r(n^)^/'^, and the 
linear biasing coefficient, (61) is evaluated at k ~ 0.1266Mpc~^. 

3.2. Power spectrum and hi 

First we estimate the linear biasing parameter according to equation (9): 

KAt)^^. (28) 



available at the website http://www-hpcc.astro.washington.edu/ 
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Fig. 5. Inverse of the linear biasing parameters of simulated halos (Table 2) and SDSS galaxies (Table 1). Crosses 
and filled circles for simulated halos (SDSS galaxies) correspond to S and L samples (blue and red), respectively. 
Dashed and solid lines indicate results based on the direct estimation of P(fe), while symbols in all panels use 
pP°'''(A:) or pP°'=(fc), equation (29). The quoted error bars for simulated halos are computed from three differ- 
ent realizations. We simply use the error bars for wedge subsamples just for reference in the case of SDSS galaxies. 

where the subscript i denotes the different subsamples of simulated halos and SDSS galaxies. 

In real-space (r-cubc in Table 2), we first construct density field of halos using the cloud-in-cell pixelization, Fourier 
transform it, angularly average Pi(k.) over the direction of k, and finally remove the shot noise contribution, l/ri/j to 
obtain Pi{k). We choose logarithmically equal bins for fc, Ailogi^k) = 1/6. The mean values {rih) of the number density 
of halos in each realization {uh) are listed in Table 2. The wavenumbers corresponding to the mean halo separation 
2tt/ {uh)^^^ roughly provide the limit of the reliability of estimating Pi{k). The dark matter power spectrum P{k) is 
estimated in two different methods. One is a direct estimate using the dark matter particles, P''™'''(fc). The other 
employs an analytical prescription by Peacock & Dodds (1996), P^^'^{k). Upper-left panel of figure 5 shows these 
results. Symbols and curves represent the results using P^^'''{k) and P''™'''(fc). Their agreement ensures the validity 
of the Peacock-Dodds prescription. 

The wavenumber dependence of is weak for k < 0.2/iMpc~^. The increase of 1/bi for k > 0.2/iMpc~^ is 
consistent with that expected from the halo finite volume exclusion effect (e.g., Taruya et al. 2001). When the 
wavenumber becomes comparable or larger than the mean separation of halos (~ 0.5/iMpc~^; c.f. Table 2), Pi{k) does 
not represent the intrinsic clustering signal properly. Therefore the increase of 1/bi against k may not be real for 
k > 0.5/iMpc~^. In this paper we are interested in linear regimes where the k dependence of bi is negligible. Thus we 
do not consider the region k > 0.2ft,Mpc~^, where the scale independent expansion like equation (4) breaks down. 

In redshift space (s-cube in Table 2), we compute Pi(fc) using the center-of-mass peculiar velocity of each halo 
adopting the distant observer approximation. The dark matter power spectrum in redshift space is computed with 
two different methods. One is a direct estimate using the dark matter particles, P^^™'^{k). The other, P^^'^{k), is 
an analytical prediction combining P^^'''(fc) and redshift distortion effects empirically. For the latter we consider the 
Kaiser effect (Kaiser 1987) and the finger-of-God effect assuming exponential peculiar velocity distribution (Cole et 
al. 1994; Cole et al. 1995; Kang et al. 2002). In this case the redshift space power spectrum is given as 
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P™''(fc), (29) 



^('^m) ~ — arctan(K,„), (30) 



B{K„,)^^[l~AiK,^)], (31) 



C(«m) = ifl-S(A.™)Y (32) 

where (Jm denotes the value of equation (27) for dark matter particles, k„i = kap/{^/2Ho), and the rms pairwise 
velocity dispersion of dark matter particles, cjp (= \/2av ~ 520km/s, cr„ is the one-dimensional velocity dispersion), is 
computed directly from N-body simulations. Upper-right panel in figure 5 displays the results for s-cube. Again bi 
is fairly scale independent in redshift space. Although the finger-of-God effect partially compensates the wavenunibcr 
dependence in r-cubc for k > 0.2ft,Mpc^^; 6i(fc) there is not reliable. The results of r-cubc and s-cubc imply that bi{k) 
is not sensitive to the redshift distortion, i.e., « b^ for k < 0.2/iMpc~"'^. In fact, the agreement between P^"'^{k) and 
P^^{k) suggests that the above feature can be explained using equation (29): 



bl^ ^ Ain„,) + lP„,B{n„,) + lPlCM' ^ ' 

where subscripts i and m refer to different subsamples and dark matter particles, respectively, and /3, indicates the 
value of equation (27) for each subsample, and Ki = kap^i/{V2Ho). The rms pairwise velocity dispersion of halos, ap^i, 
is calculated from the simulated halo subsamples through ap = \/2av; 412 km/s for S and 397 km/s for L. Equation 
(33) yields 1.04 (0.97) at A: = 0.05/iMpc~^ and 1.08 (1.01) at /c = 0.2/iMpc"^ for the S (L) subsample. These values 
explain the behavior in figure 5 and suggest that 6f and bi in linear regime agree within 10% level. 

For wedge subsamples, we determine the positions of halos in redshift space properly using the line-of-sight velocity 
component unlike s-cubc where we employ the distant observer approximation. In estimating Pi(fc), wc follow Feldman, 
Kaiser & Peacock (1994), Matarrese et al. (1997) and Verde et al. (2002), and define the field: 

F,{r) = Xw{r)[n,{r) ~ an,{r)l (34) 

where w{r) is the weight, A is a constant to be determined, nr(r) is the number density of random particles, n,;(r) is 
the number density of each subsamples, and a is the ratio of particle numbers of actual and random catalogs. In this 

— 1/2 

paper, w{r) is unity inside the survey volume, and zero otherwise. If we set A — I22 , where 

Ijk^ [ A\wHv)n^,{v) (35) 



with the mean number density for each subsample (Matarrese et al. 1997), then the power spectrum is 

(|P,(k)p)=P,(fc) + :^(l + «) (36) 

(Verde et al. 2002). The dark matter power spectrum for the wedge subsamples is calculated by the two methods 
again: the first is an analytical calculation using equation (29) and the second is a direct estimation based on equation 
(36). The results are plotted in the lower-left panel of figure 5. Symbols correspond to the estimate using equation 
(28) with P{k) evaluated from the Peacock-Dodds prescription, i.e., P^°(fc). Sohd and dashed lines use the direct 
estimate for P{k), instead, i.e., P''™(fc). The symbols and the lines deviate for k < O.l/iMpc"^ where the effect of 
complicated survey volume shape cannot be ignored. For k > 0.1/iMpc~^ the estimate of bi in wedge configuration 
well reproduces that in s-cube, upper-right panel of figure 5. 

Finally bi{k) for SDSS galaxies is computed similarly as wedge subsamples. Again wc assume the same cosmological 
parameters used in the N-body simulations. As in the wedge case, we use P{k) = P^^(fc) and P^™(fc) for symbols and 
lines, respectively. The lower-right panel of figure 5 is the result. The red (blue) sample is almost the same as L (S) in 
the lower-left panel. Lower-left panel of figure 5 is indeed in good agreement with the upper panels for k < 0.2/iMpc~^, 
which we may interpret as an indication that the SDSS results in redshift space are directly related to their real-space 
property. Wc find very similar deviations between lines and symbols at fc < 0.1, which ensures our interpretation due 
to the effect of complicated survey boundary shape. 

3.3. Bispectrum and Q 



Finally wc arc in a position to consider the three-point statistics. We calculate the bispcctra Pi(fci,fc2,fc3) using the 
same methods as Pi{k) for r-cubc and s-cube. For wedge and SDSS subsamples wc use the formula 
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Fig. 6. Qhalo fof simulated halos (Table 2) and Qgalaxy for SDSS galaxies (Table 1). Crosses and filled circles for sim- 
ulated halos (SDSS galaxies) correspond to S and L samples (blue and red), respectively. For comparison, we also plot 
the value of Qm calculated from N-body simulations (filled triangles) and from perturbation theory (PT, dashed lines) in 
r-cube, s-cube and wedge. The quoted error bars for simulated halos and dark matters are computed from three dif- 
ferent realizations. We simply use the error bars for wedge subsamplos just for reference in the case of SDSS galaxies. 



(i^,(ki)^^,(k2)^^.(k3)) = 4l \ S,:(ki,k2,k3) + ^ [P,(ki) + P,{k2) + P,{ks)] + (1 - a^)^ 

/22^ I ^33 ^33 

(Vcrdc ct al. 2002). Figure 6 plots Qi{k), the reduced amplitude of bispectrum for equilateral triangles 
B^{k,k,k) 



Q^{k) 



3Pi{k) 



(37) 



(38) 



which should be compared with figure 5. The difference of &i between the two different subsamples in each panel of 
figure 5 does not show up in figure 6. This is fully consistent with the finding of Kayo et al. (2004) for three-point 
correlation functions. Furthermore comparison among simulated halo subsamples indicates that this feature is not a 
simple outcome of the redshift distortion effect, but reflects the intrinsic correlation between bi and 62/^1 in real space 
as we discussed in section 2. 

In order to proceed further, we attempt to estimate b2{k)/bi{k) combining Qi{k) and Qm{k). We directly compute 
Qm{k) from simulation particles for r-cube, s-cube and wedge (filled triangles in fig. 6). For comparison, tree-level 
perturbation theory predicts that 

,^^2(kl,k2)P(fcl)P(fc2)+^^2(k2,k3)P(fc2)P(fc3)+^^2(k3,ki)P(fc3)P(fcl) 



g(fcl,fc2,fc3)=2- 



P{kl)P{k2) + P(fe)P(fc3) + P{k3)P{ki) 

where we follow Jain & Bcrtschinger (1994) and define i^2(ki,k2) as 
5 ^ 1 ki_k2 /^^ ^ ^ 2 (ki • k2)2 
7 2 fcifc2 \k2 ki 



F2(kl,k2) 



(39) 



(40) 
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Fig. 7. 62/61 for simulated halos (Table 2) and SDSS galaxies (Table 1). Crosses and filled circles for simulated halos (SDSS 
galaxies) correspond to S and L samples (blue and red), respectively. The quoted error bars for simulated halos arc computed from 
three different realizations. We simply use the error bars for wedge subsamples just for reference in the case of SDSS galaxies. 

Strictly speaking, this expression of i^2(ki,k2) is valid only for the Einstein-de Sitter universe, but its dependence on 
cosmology is proved to be very weak (Matsubara 1995; Scoccimarro et al. 1998). For equilateral triangles, equation 
(39) reduces to a constant value, 4/7, which is plotted as a dashed line (PT) in the upper-left panel of figure 6 for 
real space. In redshift space, one obtains Q = 0.533 by setting 7 = and (3 = Pm in the perturbation expression, 
equation (26), which is also plotted as a dashed line (PT) in the upper-right and lower-left panels. For k < 0.2/iMpc~^, 
the direct estimates (filled triangles; DM) agree with the perturbation values. In directly evaluating Q from N-body 
simulations, we use equation (8) for Qm, and a similar expression for Qi. 
Figure 7 plots &2/&1 computed from 



{b2/hUk) = bi^,{k)Q,{k)-Q„^{k). 



(41) 



For SDSS subsamples we use the value evaluated from wedge for Q,n{k). Strictly speaking, equation (41) ignores 
redshift-space distortion effects. In reality, however, we made sure that almost the same values are derived even when 
we use equation (26) instead. Therefore we use equation (41) both in real- and redshift-space. Comparison between 
figure 5 and figure 7 indicates that subsamples with larger bi tend to have larger 62/^1 ■ This is consistent qualitatively 
with analytical results shown in figure 2. 

To compare more quantitatively simulation results with analytic models, we average &i,i(fc) and (&2/&i)j(fc) over 
the range of 0.08/iMpc~^ < k < 0.2hMpc~^ . The results are plotted in the left panel of figure 8. The three curves 
show the analytic models plotted in figure 2 at z — 0, while squares correspond to the mass-averaged values for 
ellipsoidal halo model again in figure 2. Triangles, crosses, and circles, represent r-cube, s-cube, and wedge subsamples, 
respectively. If the analytic ellipsoidal halo model is exact, triangles and squares should agree within the error bars 
(we estimate from variance for three different realizations). The small differences between them may be ascribed to 
i) the inaccuracy of the higher order biasing coefficients in the halo model. The parameters {p and q) in equation 
(14) are empirically determined so as to reproduce the mass function but they do not guarantee the reliability for 
&2- ii) Weak scale dependence which is seen in figures 5 and 7. This is not expected in the analytic model, implying 
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Fig. 8. Correlation between bi and 62/^1- l^ft- analytic models and simulated halos (Table 2). Lines are analyti- 
cal predictions for spherical (solid), ellipsoidal (dashed) halo models, and peak model (dotted) as in figure 2 at z = 0. 
Squares are the mass averaged values corresponding to S and L subsamples for the ellipsoidal halo model. Other sym- 
bols are for halo subsamples. The quoted error bars are computed from three different realizations. right: SDSS 
color selected samples (Table 1; crosses), SDSS luminosity threshold samples (Table 3; triangles), HOD mock galax- 
ies (Table 3; squares), and HOD analytic prediction (circles). Note that we do not distinguish bn and B„ here. 



Table 3. HOD parameters corresponding to the SDSS luminosity threshold galaxies 



Mr^" ^min 2max N^^^^'' logip M„,in logip Afl a (iV^f^ 

-22.0 0.02 0.22 7704 13.91 14.92 1.43 1417 

-21.5 0.02 0.19 24711 13.27 14.60 1.94 7235 

-21.0 0.02 0.15 41969 12.72 14.09 1.39 28151 

-20.5 0.02 0.13 69217 12.30 13.67 1.21 75292 

-20.0 0.02 0.10 63770 12.01 13.42 1.16 141424 

-19.5 0.02 0.08 56384 11.76 13.15 1.13 246421 

-19.0 0.02 0.06 30532 11.59 12.94 1.08 365719 

-18.5 0.02 0.05 24636 11.44 12.77 1.01 514269 

-18.0 0.02 0.04 16123 11.27 12.57 0.92 745340 



its practical limitation, iii) Inaccuracy of the estimators of bi and and iv) stochasticity of biasing (Dekel 

& Lahav 1999; Taruya & Suto 2000; Taruya et al. 2001; Yoshikawa et al. 2001). The present halo model assumes 
deterministic biasing, which may affect the mean values of biasing coefficients, especially in the higher order. Given 
those realistic issues, we would interpret that triangles and squares are in reasonable agreement. 

The right panel of figure 8 plots 62/^1 against bi for SDSS galaxy subsamples, which should be compared with the 
left panel for simulated halos. In addition to the color selected subsamples in Table 1, we construct nine luminosity 
threshold subsamples (Table 3) following Zehavi et al. (2005). We construct realistic galaxy mock samples in order 
to see if the observational data can be understood from simple theoretical modeling. We employ halo occupation 
distribution approach (HOD) which assigns galaxies within simulated halos. According to the simplest version of 
HOD, mean number of galaxies within halos of mass M is set as 

(N{M)) = I ^ + (mt) ^'^ ^ (42) 
[ otherwise. 

In the above expression, the first term represents a central galaxy while the second corresponds to satellite galaxies. 
The three parameters, Mmim Mi, and a, are determined to reproduce the observed sample of galaxies. In practice, 
we adopt the values (Table 3) fitted by Zehavi et al. (2005) for the SDSS luminosity galaxies (see their Table 2). 
Their fit assumes the density profile proposed by Navarro, Frenk & White (1996) for galaxy distribution in each halo, 
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and satellite galaxies within each halo are assigned following the Poisson distribution of the mean value of equation 
(42). We construct mock galaxy samples from our simulated halo catalogs using the routine of Skibba et al. (2006) 
which simultaneously takes care of the redshift-space distortion effect (see also Kang et al. 2005). In Table 3, N^^^^ 

and {N^°'^^) denote the numbers of SDSS galaxy samples and the mean numbers of mock galaxy samples. Since the 
galaxy samples of Zehavi et al. (2005) are based on different galaxy catalog (SDSS DR2; Abazajian et al. 2004) from 
this work, the galaxy numbers of our luminosity threshold samples are about twice larger than theirs. We simply 
distribute mock galaxies in the simulation cubes and do not consider the survey boundary effect, so the numbers of 
mock galaxy samples do not correspond to those of SDSS galaxy samples. 

According to the current HOD approach, we evaluate the biasing coefficients in two ways. One is an analytical 
estimate (HOD analytic prediction) which is based on equation (25): 



where we use the Sheth-Tormen mass function, equation (14), for nhaio(^^, ^ = 0), and equation (20) for bn{M,z = 0), 
and Mmax is set as the maximum mass of our simulated halos. The other is the direct evaluation of the biasing 
coefficients from the HOD mock galaxy samples using equations (28) and (41). The right panel of figure 8 plot these 
results for SDSS color selected samples (crosses), for SDSS luminosity threshold samples (triangles), for HOD analytic 
prediction (circles) and finally for HOD mock galaxies (squares with error bars). We note a few interesting features 
in the plot, i) The similarity between the color selected and luminosity threshold samples. While they are based on 
different selection criteria, the resulting &2/^i ~ bi correlation seems to be roughly the same, ii) The HOD analytic 
prediction reproduces the values of 62/^1 and bi for each luminosity threshold galaxy samples. This is surprising 
since the fitting procedure is designed to reproduce their two-point correlation functions alone, iii) The discrepancy 
between the HOD mock and analytical results may come from the difference between Shcth-Tormcn model and our 
halo samples in rihaio and/or 62. We made sure that nhaio of our halo samples is in good agreement with equation (14). 
On the other hand as the left panel of figure 8 indicates, 62 for our simulated halos is systematically lower than the 
analytic prediction, equation (25). Thus the agreement between SDSS data and HOD models needs to be interpreted 
with caution. We also compute the values of &2/&1 and hi of HOD mock galaxy samples in real space. We, however, 
do not plot the results, since they are almost the same as those in redshift space. 

Even with the above subtlety we find a fairly generic correlation between 62/&1 and &i, which is the most important 
result of the current study. We note here that Gaztafiaga et al. (2005) performed a related analysis using 2dF galaxy 
data. They found that bi = 0.93]J^p pg and 62/&1 = — 0.34^'^p Qg . These values are in good agreement with the right 
panel of figure 8. Incidentally they speculated a crude correlation of 62/^1 ^ ^ 1-2, which is indeed consistent with 
our finding on the basis of systematic results from analytic, and numerical simulations and SDSS galaxy data. 

4. Summary and Discussion 

We have found a fairly generic correlation between linear and quadratic biasing coefficients, 61 and 62 /&i, using 
a variety of different methodologies; perturbative expansions of peak and halo biasing models, N-body simulations 
of halos, SDSS galaxy data analysis, and the corresponding halo occupation distribution predictions (analytic and 
mock). The presence of such correlations was suggested earlier by the previous finding that the normalized three- 
point correlation functions of SDSS galaxies in redshift space follow the hierarchical relation approximately, Q = 0.5 ^ 
1.0, despite the robust morphological, color and luminosity dependences of the corresponding two-point correlation 
functions (Kayo et al. 2004). We have derived the bi - 62 correlation explicitly and showed for the first time that it 
indeed explains the observed behavior of Q for equilateral triangles calculated from SDSS galaxies in linear regimes. 

The major findings of the present paper arc summarized as follows. 

• Even with the presence of redshift distortion and complicated survey shape effects, &i can be accurately estimated 
in fc < 0.2ft,Mpc~^. So 61 estimated from SDSS galaxies is expected to reflect the true value in real space. 

• The values of Qb for equilateral triangles from simulated halos, mock galaxies based on HOD model, and SDSS 
galaxies do not show any clear dependence on 61. The independence on bi is not an artifact from the redshift- 
space distortion contamination, but is a consequence of the intrinsic correlation between 62/^1 and bi. 

• The HOD model, equation (42), whose parameters are determined so as to reproduce the observed two-point 
statistics alone, seems to be also successful in predicting the value of 62. 

There are a few remaining tasks which should be done following the present result. First, it is interesting to compare 
the current analysis in Fourier space with that in configuration space. As long as nonlinear effect is negligible, Qm 
and Qb in k-space are expected to be the same as those in configuration space. Nevertheless we would like to make 




0) 



(n = l,2). 



(43) 
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sure of it, and to further explore the connection between bispectra and the three-point correlation functions of halos 
and galaxies. Second, the success of HOD to reproduce the correlation between 62/61 and 61 at the current level is 
very encouraging. Naturally it is important to further improve the HOD model. Third, the present paper considers 
only equilateral configuration of Fourier space triangles for Q. In reality, however, the different shape of triangles has 
a wealth of information. While the 61 dependence of Q barely vanishes for equilateral triangles due to the correlation 
between 61 and 62/61, this is not the case for arbitrary triangle shapes. Therefore it is interesting to explore if the 61 
dependence of Q for the other triangle shapes and to see if the results are consistent with the present analysis and/or 
to constrain higher order biasing coefficients more tightly. The work along the line is now in progress. 
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